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Abstract 

The predicted value of the higgs mass run is analyzed assuming 
the existence of the fourth generation of leptons (N, E) and quarks 
(U, D). The steep and flat directions are found in the five-dimensional 
parameter space: mjj, mu, rno, m^, the- The LEPTOP fit of the 
precision electroweak data is compatible (in particular) with run ~ 
300 GeV, m N ~ 50 GeV, m E ~ 100 GeV, m v + m D ~ 500 GeV, and 
\mjj — mo\ ~ 75 GeV. The quality of fits drastically improves when 
the data on b- and c-quark asymmetries and new NuTeV data on deep 
inelastic scattering are ignored. 

It is well known that in the framework of Standard Model the fit of 
electroweak precision data results in prediction of light higgs, the central 
value of its mass being lower than the direct lower limit set by LEP II |l[ . One 
possible way to raise the predicted value of m# is to assume the existence of 
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fourth generation of leptons and quarks, [0, [| . Implications of extra quark- 
lepton generations for precision data were studied in a number of papers 
Leptons of fourth generation (E,N) should be very weakly mixed 
with the ordinary ones, while in quark sector (U,D) mixing is limited only by 
unitarity of 3 x 3 CKM matrix. In particular it was noticed in ref. || that the 
predicted mass of the higgs could be as high as 500 GeV. That conclusion was 
based on a sample of 10.000 random inputs of masses of fourth generation 
leptons and quarks. However the sets of the lepton and quark masses were 
presented independently (see Fig. 7 in ref. Thus it is not clear how they 
were combined. 

In this letter we try to develop a systematic approach to the problem 
by using our LEPTOP code || to find steep and flat directions in the five- 
dimensional parameter space: m#, mu, itld, the, tun- For each point in 
this space we perform three-parameter fit (rrit, a s , a) and calculate the \ 2 °f 
the fit. It turns out that the xLin depends weakly on mu + mo and m#, 
while its dependence on mjj — mp, m E and is strong. We limit ourselves 
to the values of larger than 50 GeV because according to experimental 
data from LEP II on the emission of initial state bremsstrahlung photons, 
m N > 50 GeV at 95% c.l. [0, [OJ. 

We analyzed Summer 2001 precision data (ref. [jlj which are also given 
in the Table 1 in ref. ||). Figures 1-4 show xLin (crosses) and constant \ 2 
lines corresponding to A% 2 = 1,4,9, 16, ... on the plane m N ,mu — m D for 
fixed values of mu + m D = 500 GeV, m H = 120 (Figs. 1 and 3) and 500 GeV 
(Figs. 2 and 4) and m E = 100 (Figs. 1 and 2) and 300 GeV (Figs 3 and 4). 
We also performed fits for mn = 300 GeV. 

The above choice of masses is based on a large number of fits covering 
a broad space of parameters: 300 GeV < m v + m_o < 800 GeV; GeV 
<mu-m D < 400 GeV; 100 GeV < m E < 500 GeV; 50 GeV < m N < 500 
GeV; 120 GeV < m H < 500 GeV. Concerning quarks, mu + m D is bounded 
from below by direct searches limit, while from above by triviality arguments. 
Since x 2 dependence on mu + m E is very weak, our choice of intermediate 
value mu + rn D = 500 GeV represents a typical, almost general case. For 
this choice \mu — mo\ can not be larger than ~ 200 GeV because of the 
mentioned above direct searches bound. 

Concerning charged lepton, its mass is taken above LEP II bound. We 
present fits at two values of (100 GeV and 300 GeV) and one can see 
how fit is worsening with m E going up. 

Concerning the value of mn, we vary it from the lower LEP II limit up to 
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triviality bound and since the dependence of observables on m# is flat, one 
can get x 2 behaviour from two limiting points: m# =120 and 500 GeV. 
For rriE = 100 GeV we have the minimum of x 2 & t m N — 50 GeV and: 



for 


m H 


= 120 GeV: 


\m v 


- m D 


~ 50 GeV, 


Xmm/ n d.o.f. 


= 20.6/12 


for 




= 300 GeV: 


\m v 


- m D 


~ 75 GeV, 


Xvam/ n d.o.f. 


= 20.8/12 


for 


m H 


= 500 GeV: 


\™>U 


- m D 


~ 85 GeV, 


Xnan/ n d.a.f. 


= 21.4/12 



Thus we have two lines (mu > and mjj < mo) in the (mu — mp, m#) 
space that correspond to the best fit of data. Along these lines the quality 
of the fit is only slightly better for light higgs (m# ~ 120 GeV) than for the 
heavy one (m H ~ 300 - 500 GeV). 

Note that the n d _ _f is 12, unlike the case of the Standard Model, where 
it was 13 (ref. ||). This change occurs because in the present paper m# 
is not a fitted, but a fixed parameter (hence 13 becomes 14), while 
and mjj — mn are two additional fitted parameters (hence 14 becomes 12). 
(As is well known, rid.o.f. is equal to the number of experimentally measured 
observables minus the number of fitted parameters.) 

For rriE = 300 GeV we have the minimum of x 2 & t my ~ m D — 25 GeV 
and: 



for m H = 120 GeV 
for m H = 300 GeV 
for m H = 500 GeV 



m N ~ 200 GeV, xL/% ./. = 23.0/12 
m N ~ 170 GeV, xljn d . . f . = 24.0/12 
m N ~ 150 GeV, xLJn d . . f . = 24.4/12 



Thus, the best fit of the data corresponds to the light tre — 100 GeV 
and Tn 7y — 50 GeV. The significance of light (around 50 GeV) was 
first stressed in p. Increase of leads to the increase of m N and to fast 
worsening of Xmin- 

Although inclusion of one extra generation improves the quality of the fit 
(compare x 2 /n d .o.f. = 23.8/13 for the SM from [|] and X 2 mi Jn d .o.f. = 20.6/12 
from Fig. 1) it remains pretty poor. The poor quality of the fit is due to 3.3cr 
discrepancy in sf = sin 2 6^ extracted from leptonic decays and from A b F c B 
I T|. If one multiplies experimental errors of A b FB and A C FB by a factor 10, one 



gets good quality of SM fit JTT], g| but with extremely light higgs, having only 
a small (few percent) likelihood to be consistent with the lower limit from 
direct searches. We prove that the fourth generation allows to have higgs as 
heavy as 500 GeV with a perfect quality of the fit: Xmin/ n d.o.f. = 13/12, if 
one uses old NuTeV data (see caption of Fig. 2). 
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To qualitatively understand the dependence of mjj — mj) on m# in the 
case of m E = 100 GeV at xLm ^ us reca U how radiative corrections to the 
ratio mw/mz and to qa and R = gy / ' Qa (the axial and the ratio of vector and 
axial couplings of Z-boson to charged leptons) depend on these quantities 
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11 2 \ / 16 G 2 m v ~m D 

2,1 / 3 m z \ 2 m u -m D 

' 9 / \ 9 my+mo 



(1) 



where i — m, A, R, while s 2 ~ 0.23. Corrections to other observables can be 
calculated in terms of 5V % . In the vicinity of Xmm t ne third term in brackets 
is much smaller than the second one. Hence the smallness of the left-right 
asymmetry of the plots of Figs. 1,2. Since ^-s 2 ~ s 2 + | ~ s 2 , the increase 
of rrin is compensated by increase of \mu — m D \ and we have a valley of Xmm- 
Captions of Figs. 1 and 2 reflect recent change in NuTeV data (from 



m w = 80.26 ± 0.11 GeV to m w = 80.14 ± 0.08 GeV 0) which results 
in drastic worsening of the fit even in the presence of the fourth generation. 

Thus we see that the 4th family scenario is better than the Standard 
Model, because the latter can produce good fit only when the mass of the 
higgs is much lower than the lower limit of LEP II, even when experimental 
data on heavy quark asymmetries and new NuTeV data are ignored. 



Note that originally introduced in \Ti\ parameters S,T,U are not ade- 
quate for the above analysis, because they assume that all particles of the 
fourth generation are much heavier than mz, while in our case the best fit 
corresponds to ~ mz/2. In the paper modified definitions of S and U 
were used in order to deal with new particles with masses comparable to m z . 
However, let us stress that both original and modified definitions of S, T 
and U take into account radiative corrections from the "light" 4th neutrino 
only approximately, while the threshold effects, that are so important for 
m N ~ 50 GeV, can be adequately described in the framework of functions 
V\ 

In conclusion let us stress that in the framework of SUSY with three 
generations radiative corrections due to loops with superpartners also shift 
upward the mass of the higgs in the case of not too heavy squarks (300-400 
GeV, see Table 1 in UT5| ) or light sneutrinos (55-80 GeV, see fL6"f ). 
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Figure 1: Exclusion plot on the plane mN,mu — ran for fixed values of 
itlh = 120 GeV, my + mo = 500 GeV and itle = 100 GeV. Xmm shown 
by two crosses corresponds to x 2 / n d.o.f. = 20.6/12. (The left-hand cross 
is slightly below tun = 50 GeV.) Borders of regions show domains allowed 
at the level Ax 2 = 1,4,9,16, etc. The plot was based on the old NuTeV 
data. The new NuTeV data preserve the pattern of the plot, but lead to 
Xmm/ n d.o.f. — 27.7/12. If A b FB and A C FB uncertainties are multiplied by factor 
10 we get X 2 m Jn d .o.f. = 19.1/12 for new NuTeV, and X 2 m Jn d . oJ . = 11.3/12 
for old NuTeV with practically the same pattern of the plot. 
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Figure 2: Exclusion plot on the plane m^^mu — vtld for fixed values of 
m H = 500 GeV, my + m D = 500 GeV and uie = 100 GeV. Xmin shown 
by two crosses corresponds to x 2 / n d.o.f. = 21.4/12. (The left-hand cross 
is slightly below = 50 GeV.) Borders of regions show domains allowed 
at the level A^ 2 = 1,4,9,16, etc. The plot was based on the old NuTeV 
data. The new NuTeV data preserve the pattern of the plot, but lead to 
Xmm/ n d.o.f. = 28.3/12. If A b FB and A C FB uncertainties are multiplied by a 
factor 10, we get X 2 m iJn d .o.f. = 21.2/12 for new NuTeV, and xLiJ 'n d .o.f. = 
13/12 for old NuTeV with practically the same pattern of the plot. 



8 




-200 -150 -100 -50 50 100 150 200 



niy-mp (GeV) 



Figure 3: Exclusion plot on the plane m^^mu — m D for fixed values of 
m H = 120 GeV, my + m D = 500 GeV and m E = 300 GeV. Xmm shown by 
cross corresponds to x 2 / n d.oj. = 23.0/12. Borders of regions show domains 
allowed at the level A% 2 = 1, 4, 9, 16, etc. 
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Figure 4: Exclusion plot on the plane m^^mu — m D for fixed values of 
m H = 500 GeV, my + m D = 500 GeV and m E = 300 GeV. Xmm shown by 
cross corresponds to x 2 / n d.oj. = 24.4/12. Borders of regions show domains 
allowed at the level A% 2 = 1, 4, 9, 16, etc. 
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